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FOREWORD 


' The free electron laser instability for a solid relativistic electron beam 
propagating in combined transverse helical wiggler and uniform axial guide 
fields is investigated within the framework of the linearized Vlasov-Maxwell 

equations. Ht is assumed that v/y b <<1, where v is Budker's parameter and 

2 -'v 

Y b mc is the electron energy. 'Stability properties are investigated for the 
choice of equalibriurn distribution function in which all electrons have the 
same value of the linear combination of transverse and helical invariants./ 

C 1 - 2y b mu) b C h = const., and a Lorentzian distribution in the axial invariant C z> 
(Here w b is a constant.) The instability growth rate is calculated including a 
determination of the optimum value of the ratio of beam radius to conducting 
wall radius (Rq/R c ) for maximum growth. It is found that the maximum growth 
rate for a solid electron beam is comparable to that for a hollow beam with 
similar parameters. Moreover, the introduction of a small axial momentum spread 
(^/Y b irar^ff few percent) significantly reduces the instability growth rate, pr- 
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1. INTRODUCTION 

1-9 

In recent years, the free electron laser instability has been 
extensively Investigated with particular emphasis on applications 
to intense microwave generation. For the most part, previous theoretical 

analyses of this instability have been carried out either for uniform 
3-6 

density beams with infinite transverse dimension, or for annular 
7 8 

electron beams. ’ The present paper examines the influence of finite 
radial geometry on the free electron laser instability for a solid 
electron beam propagating in combined transverse helical wlggler and 
uniform axial guide fields. The analysis is carried out within the 
framework of the linearized Vlasov-Haxvell equations, including a 
determination of the optimum value of beam radius R^ for maximum 
growth rate. 

The present analysis is carried out for an infinitely long relativ¬ 
istic electron beam propagating in the combined transverse wiggler and 
uniform axial guide fields described by Eq. (1). Equilibrium and stability 
properties are calculated for the specific choice of electron distribution 
function (Eq. (5)], 

f b ' i(Ci " 2 Yb ra "b C h “ 2 Y b ' n T 1 )G(C z ) , 

where n^, y^. and T x are positive constants, C i( C^, and C z are 

the transverse, helical, and axial invariants*^ defined in Eqs. 

(6) - (8), and the axial distribution function is normalized according 
to dC z G(C z ) ■ 1. Equilibrium properties are investigated in Sec. II, 

and stability properties are examined in Secs. Ill and IV, assuming 

2 

that v/y^ « 1, where v is Budker's parameter, y^mc is the characteristic 

electron energy, and m is the electron rest mass. j 

I 
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In Sec. Ill, making use of the linearized Vlasov-Maxwell equations, 
we obtain the coupled eigenvalue equations (33), (37), (40), (45), and 
(46) that describe free electron laser stability properties in 
circumstances where the perturbed transverse fields can be approximated 
by the vacuum waveguide fields. For short wavelength perturbations, 
the axial component of the perturbed longitudinal field can be approximated 
by (Eq. (49)], 


£ X<'>' 


\ 0 , ot 


r/R 0 ) , 


0 < r < R„ 


otherwise 


where Jj(x) is the Bessel function of the first kind of order l, 

6. , is the s'th zero of J (B, ,) - 0 and $ , is a constant. In 

Sec. IV, substituting Eq. (49) into the coupled eigenvalue equations, 
we obtain closed algebraic dispersion relations for the transverse 
electric (TE) and transverse magnetic (TM) polarizations. 


Introducing the normalized dimensionless function [Eq. (59)], 

2 J £ (x) 

G ts- (x) ■ S,.- tt~ r—2 * 

(X 6 l, s’> 

it is shown in Sec. IV that the coupling between the longitudinal and 


transverse perturbations Is proportional to G ,(o ., E n /R ) 

£S fc'i|8 U C 

for the TE mode, and to G. R n /R ) for the TM mode. Here 

£s £+l,s 0 c 

“t+l.S and e £+l,s are the sthzeroes ° f j ;+lS+l,s ) ■ ° and J i+l (B £+l,s ) “ °* 
respectively, R c is the radius of the outer conducting wall, and the prime 

denotes (d/dx)J (x). Assuming that the maximum of the function C , (x) 

£ £S 

occurs at i ■ x ,, we note that the maximum instability growth rate 
Is 

occurs at a value of R„/R given by R n /R » x. for the TE 

0 c 1 0 c is l+l,s 

mode, and by Rq/R c ” x^ g ,/8 t+ ^ s Bor t * le ® mode. This result is 

g 

different from that obtained for a hollow electron beam. 


8 
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A detailed numerical analysis of the TE mode [Eq. (62)] and TM 

mode [Eq. (63)] dispersion relations is presented in Sec. IV. Two 

features are noteworthy from the numerical analysis. First, for the optimized 

value of R n /R , the instability growth rate for the TM mode is comparable 
u c 

to that for the TE mode. Moreover, the growth rate is reduced substantially 
by introducing a small amount of axial momentum spread (A/y^mc £ 0.01). 


9/10 
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II. EQUILIBRIUM THEORY AND BASIC ASSUMPTIONS 


The equilibrium configuration consists of a relativistic electron 
beam propagating in the combined transverse wiggler and uniform axial 
guide fields described by 


B = -SBcos(0 - k 0 x)£ r + 6Bsin(8 - k 0 z)e e + B^ 


where Bq and 6B are constants, and kg is the axial wavenumber of the 
helical wiggler field. In Eq. (1), cylindrical polar coordinates (r,6,z) 
are used, with z-axis along the propagation direction, and and 

are unit vectors in the r-, 8-, and z-directions, respectively. In 
the present analysis, we assume that the axial wavenumber of the helical 
wiggler field is sufficiently large that 

"l-- i- ~\ « U 0 H , (2) 

l w 0 “ V V 00 

where wg " k^V^, “c = e ®0^ Y b mc * S t * ie e ^ ectron cyclotron frequency, 

2 

Rq is the characteristic beam radius, Y^mc is the characteristic 

2 1/2 

electron energy, c is the speed of light in vacu o, » c(y fa - 1) /y^ 

is the mean axial velocity of the electron beam, and -e and m are the 
electron charge and rest mass, respectively. 

It is also assumed that 


v/yv « 1. 


where v *> N,e /me is Budker's parameter, 


, 2 * ,R_ 

d6 i 


dr r n b (r, 6 - k Q z) , 


is the number of electrons per unit axial length of the beam, n^Cr, e_ k^z) 
is the equilibrium electron density, and R c is the radius of the conducting 
wall. The inequality in Lq. (3) indicates that the beam is very tenuous, 


. . . »V : —o in — ii<i a r& KTmKn * 


11 
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In the parameter regimes of practical interest for free electron 

laser applications, the axial distribution function G(C^) is strongly 

peaked about « Y b mV b " const, » with characteristic half-width 

AC « y mV . Moreover, in the present analysis. It is also assumed 
Z D D 

that the axial motion is nonresonant with 


.2 2 , 2 

k 0 v x * “c ’ 


where v z * p z /ym is t ^ le axial velocity of a typical beam electron. 


( 10 ) 


We therefore approximate Eq. (8) by 


in 6B 


10 


C z “ P z + (l -J TbI [p r cos(e - k o*> ■ P 6 sin(e - k o z)] - 

0 c 0 


( 11 ) 


Making use of Eqs. (6), (7), and (11), it is straightforward 


to show that the combinat ii>o C A • In Kq. (5) ran be expressed as 


10 


Ci - 2 Ybin c h « (p r - ) cos(e - k 0 *)J 

4 

+ ( p e - Y b m V + ^(4^) Sln(6 ' k ° z) ] 


( 12 ) 


+ y 2 b * 2 + 0 ir, 0 - k Q z) , 


where tu Q - k Q V b> “ c “ eB o^ Y b mc ’ and the effective potential ^(r, 0 - kyZ) 
is defined by 


f u c " “b\ r_ SB 
“c 1 ^ - wj k 0 B 0 

. ( 3LL%Y 4 f V? + - ( C - Y b nV, 

\“o ' “c/ J ' B o' b V\ “0 / b r b b 


^0 (r ’ 6 " k 0 z) " (u *b“c ■ %^ r + 2 

k 2 2 


sin(e - k Q z) 


(13) 


As a simple example, we consider an axial distribution function 


in which all electrons have a same value of C , i.e., 

z 
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where the mean radius Rq is defined by 


y. m V - u / 2 \ B„ / b 
b 0 c w n 0 


( “b“c ~ V 

and use has been made of Eq. (2). Additional general equilibrium 
properties associated with the distribution function in Eq. (5), 


including helical distortions of the beam equilibrium for finite 6 B/Bq, 
are discussed in Ref. 10. 


14 
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III. LINEARIZED VLASOV-MAXWELL EQUATIONS FOR A TENUOUS BEAM 


In this section, we make use of the linearized Vlasov-Maxwell 
equations to investigate the free electron laser stability properties 
of a relativistic solid electron beam described by the equilibrium 
distribution function in Eq. (5). We adopt a normal-mode approach in 
which all perturbations are assumed to vary with time and space 
according to 

6\J/(x,t) = l i^ n \r)exp{l[£6 + (k + nk Q )z - wt]) , (17) 

i,n 

where lmm > 0. Here, uj is the complex eigenfrequency f k + nk^ is the 
axial wavenumber, and l and n are integers. Moreover, it is also assumed 
that the perturbations are close to resonance with 


|w - (k + nk Q )V b | lo^, w c , (18) 

where ujg » k Q V b and u> c « eB Q /Y b mc. 

The Maxwell equations for the perturbed electric and magnetic field 
amplitudes can be expressed as 


Z » £(x> “ i(w/c)g(x) , 

V * B(x) - (4i./c)J(jc) - i (w/c)E(x) , 


(19) 


where 


” ' e 

is the perturbed current density. 

r° 

f b (x,^) * e dx exp(-iut) 
J — oo 


} d3 P Z £ b ( ?*«> 


In Eq. (20), 


E(x’) + 

% a 


v' x B(x') 



( 20 ) 


( 21 ) 


15 


j. 
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is the perturbed distribution function, i - t' - t, and the particle 
trajectories ^'(t 1 ) and jj'(t') satisfy djjc'/dt' - and djg'/dt' « -e^' * ^/c, 
with "initial" conditions jg'(t' ■ t) ■ and v‘(t' ■ t) » 

Within the context of Eqs. (3) and (18), the perturbed distribution 
function in Eq. (21) can be approximated by 


f b (^t,jg) ~ ~ ~ J dr exp(-iur)|2 ^ymiui(v' • Jg) 

-PzU’-bl)) ^2 f S + (*' * h i)\r, f °} • 


( 22 ) 


2 22 1/2 

where y » (1 + p /m c ) , and use has been made of Eq. (19). To lowest 


order, the axial motion of an electron is free-streaming with 


10 



ym 


t) . 


(23) 


Moreover, within the context of Eq. (18), on the right-hand side of 
Eq. (22) we retain contributions to and in the orbit integral of 
the form 8 


U'c j g 

v' * V -— cos(6 - k n z - k_.v t) , 

r z cu_ - a) B n 0 0 z ’ 

0 c U 


(24) 


and 


z w 0 “ “c ®0 


C sin(6 - kQ* - k o v z T ^ 


(25) 


Finally, since the oscillatory modulation of the radial and azimuthal 
orbits is small amplitude (Eq. (2)], we approximate 


= r , e' = e , 


(26) 


in the arguments of the perturbation amplitudes on the right-hand side of 
Eq. (22). 


16 
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Substituting Eqs. (23) - (26) into Eq. (22), we obtain the perturbed 
distribution function 


f b ( ^’? ) “ E fb" )e *P UU6 + (k + nk Q )z)) 


i.n 


iec l eXpll ' t6 t (k + nk ° )2l 7 X C £ (n) 
“ u ' < k nk 0 )v z \ “ 2 2,1 


(r) + 


e6B 


2ymc 2 kg “0 


(27) 

wnere the function X n ,(j>,ii,,k) is defined by 


3f h 8f b 

X n '(£"“» k ) “ 2(ym4i - (k + n'kjj)p z ] — j + (k + n'kg) -- 


(28) 


3 Pi 


and 8 * v /c. In Eq. (27), the term proportional to X is the 

z z n 


longitudinal portion of the perturbed distribution function. 

Similarly, the terns proportional to X^ +1 and X j in Eq. (27) are the 
transverse electromagnetic portions of the perturbed distribution function. 


Consistent with Eq. (18), the eigenfrequency ui can be approximated 


7 2 2 

by u = (k + nk o^ V b' We therefore approximate u /c - (k + nkg + kg) by 


ui 2 /c 2 - (k + nkQ + kg) 2 » 


(k + nk_) _ 

- ~ + 2k 0 (k + nk Q ) + k 2 


< 0 , 
(29) 


for k + nkg > 0. Evidently, Eq. (29) indicates that the n + 1 mode in 
Eq. (27) is a non-propagating wave in a vacuum waveguide. Without loss 
of generality, for a tenuous beam, we therefore assume 


£(n+l) 

r,i-l 


(r) 


p(n+l) 

8 , 1-1 


(r) 


0 . 


(30) 


in the subsequent analysis. Making use of Eq. (30), ^(x,^) in Eq. (27) 
can then be expressed as 


17 




NSWC TR 81-145 


iec exp{i[t8 + (k + nk^z]) 

u ^ u - (k + nk.)v 

1 ,n 0 2 


VB E* n J< r) 

| n z z i x 


+ AA n-l^r,l+l (r) ' i ^0,£+l (r) 1 


}• 


(31) 


(32) 


where the dimensionless parameter A is defined by 

e6B U 0 

A " , 2, u 0 - u ’ 

2y^mc k Q 0 c 

and use has been made of the approximation y - y^, which is consistent 
with Eq. (18). 

From Poisson's equation, V • £(jjO “ 4rp(x), and the 
Maxwell equation (19), we obtain the differential equation, 


[ V ' + 4 - <k + nk 0 ) 2 ]i£> 


4,i(k + nk 0 ) 

( r ) = --- (r) 


(33) 


for the axial (longitudinal) component of the perturbed elecLric 


field £<">. In Eq. (33), pf°(r) 


-e i 


d^pf^"^ is the perturbed charge 


density, V 2 i r -1 (3/3r) (r3/3 r) - i 2 /r 2 , and use has been made of J^(r) 
V b p£ n ^(r). In the tenuous beam limit (Eq. (3)], the transverse field 
components (r) in Eq. (31) can be approximated by the vacuum 

g 

waveguide fields. In this context, the present stability analysis 
utilises the vacuum transverse electric (TE) and transverse magnetic 
(XM) waveguide modes as a convenient basis to represent the general 
electromagnetic field perturbation ^"^^(r), which is determined from 


^ - (k + nk Q 


k 0 > 




(34) 


“ l 1 32 fe 2?£+l ( ^ ) 1 C Iz * ? i ®2,*+l ( ? ) 




Making use of Eqs. (19) and (34), and neglecting the perturbed 
current density, the vacuum waveguide fields can be expressed as 


18 
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l £>> - 


for the TE mode, and 
i(n-l) 


loR 


(35) 


b - J . K+i > . 


“t+l.s *+l.*V a »+l.. r/K c 


2,£+l 

Un-l) 


(r) “ £ W,s J t+ l ( ' WlS r/R c> • 


(3b) 


F —- «(«-!).. _ . (k + nk 0 - k 0> R c 

r.*+l ( 6, t+ l (r) 1 ^- L t + l,sWl,s r/R c> - 

for the TO mode. In Eqs. (35) and (36), b £+1>s and c £+ls are constants. 

Jj£* (x) is the Bessel function of first kind of order fc' and a 

’ fc+l,s 

and 6 £+1>s are the sth roots of j; +1 (a £+1>s ) - 0 and 3^(8^) - 0, 
respectively. Here the prime (') denotes J £lJ (x) « (d/dx) J £+£ (x). 

Attei auue straightiorwai .1 algebraic uianipulation of Eqs. (19), (35), and 
(36), it can be shown that 


f 2 a 2 

^-q i .i _ -,2 t+l,s 

| c 2 (k+nk 0'V 

i R 

c 


b j ( -iM» s \ 
£+l,s J £+i\ r ) 


-£{h !>>;"(„} . 


(37) 


for the TE mode, and 


^ 4 - (k + nk. - k ,) 2 - g ^ 1 »£ r J ( bl + 1 > sT ) 
c 2 0 V r 2 f X+l,s J £+lV R ) 




(38) 


’ M (k + nk o * k o ) ^+i l)(r) - 7 3 ^m (r) ) ■ 

for the TO mode. Moreover, making use of the continuity equation, 

i(k + nk n - k„)j (n_1 | = i L. [r j( n -l) l + Ml ± 1) j(n-l) 
1+1 0 k 0' j z,£+ 1 r Jr rJ r,i+l' + r J e,l+I* 

(39) 


19 
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the approximation = V^^j^fr) (consistent with Eq. (3)], 

and approximating k + nk^ = k Q /(l - V^/c) on the right-hand side of 
Eq. (38), we find that tq. (38) can be expressed as 


^2 - (k + nk Q - 
c 


k 0 ) 


" i+1,s 
2 


£ j /fi+iial 

t+i.s i+n r 


An Is _ 

rc |3r 


1 rJ 


(n-1) 
r, i+1 


(r» + i(^)J^(r)} . 


(40) 


for the TM mode. 

For convenience of notation in the subsequent analysis, we intro¬ 
duce the effective susceptibility, 


X n n >(w.k) = 4ne 2 d 3 p 


A ^,u),k) 

U) — (k + nk Q )v z 


(41) 


Moreover, to simplify the present analysis, we also assume that the 
beam rotation is slow with 


“b << V “o • (i2) 

Within the context of Eq. (42), we can show from Eq. (12) that the 
equilibrium distribution function is an even function of 


p r - 2y b mcAcos(8 - kgZ) (43) 

and 

p Q + 2Y b mcAsin(6 - k Q z) , (44) 


for the beam rotations satisfying << w c> ujq. Making use of 
Eqs. (31), (41), and (43), the perturbed charge and current densities 
are given by 


‘(n-1) 

e,i+i 


(r) 



4irw 


Gj(u,k,r) 


-(n-1 ). . 
r,i+l U ' 


icAp^ n ^(r) 


(45) 


where the function Gg(u>,k,r) is defined by 


20 
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C^w.k.r) 


* 6 E (n) 

*n, n z 2,1 


(r) 


+ A *n,n-l [i ^ (r) 


^(n-1) . . . 

XL e,f+l lr ' J 


(46) 


Equations (33), (37), and (40), when combined with Eq. (45), constitute 
one of the principal results of this paper and can be used to investigate 
stability properties for a broad range of system parameters. Moreover, 
in limiting cases, the dispersion relation for the free electron laser 
instability can be obtained in a closed form (Sec. IV). 


21/22 
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XV. FREE ELECTRON LASER STABILITY PROPERTIES 


In this section, simplified expressions are obtained for the longi¬ 
tudinal perturbations in Eq. (33),and the results are used to derive the 
dispersion relation for several values of azimuthal harmonic number t. 

The present analysis assumes short wavelength perturbations with 

|q*| = |(k + nk 0 ) 2 - u 2 /c 2 | » 1/R 2 . (47) 

Moreover, for u = (k + nk^)V^ and k + nk^ - k^/(l - V^/c), the inequality 
in Eq. (47) can be expressed in the equivalent form, 

(1 + V b /c) 2 Y ^k 2 R 2 » 1 , (48) 


which is readily satisfied in the parameter regimes of present experimental 
interest. As shown in Appendix A, for short wavelength perturbations 
satisfying Eq. (48), the axial component of the perturbed electric field 
E^ n ^(r) in Eq. (33) can be approximated by 

Z , l 


p(n) 

z,£ 


(r) 


0 * 


e/Rq) , o _< r 
otherwise . 


< K, 


(49) 


In Eq. (49), 6 , is the s'th root of J (G. ■) = 0, and it , is a 

St f S A 1 J S JL y S 

constant. 

Substituting Eqs. (35) and (49) into Eq. (46), multiplying Eqs. 

(33) and (37) by rJ £ (6 £ g ,r/R 0 ) and rJ i+l^ u £+l s r / F ' c >» respectively, 
and integrating from r * 0 to r » R^, we obtain two homogeneous equations 
relating the perturbation amplitudes , and b . For the TE 

S AT 2. y S 

mode polarization, these are 

Vf C' ' ‘"o - *> (7 A,. *\*^) j2 . (\’- ) 

b 0 (50) 
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and 


(k + nk_)R 
_ 0 c 

Y b a t+l,s 




i+l 


fR 

c 

* S J 0 


dr r 


~T (k + nk- - kj 2 - + A 2 ) 


n,n-l 


,2 I a l+ l,s r 


1+1 




c 

(51) 


irHsjiq*v.-| 0 c '•<«» - ~ s ^ I ) J t+i(- ! 5^ £ ') • 


where •(x) is the Heaviside step function defined by 


«(x) 


/ 1 • * ‘ 0 • 

[ 0 , otherwise 


(52) 


Similarly, for the TM mode polarization, we obtain 

. 2 


and 


r R 

•*!; 


(53) 


♦.,.■1;^« *<»o - «> (t «... * t ) 

^ »» fc.11 . « 1 itr . 

y*’ ( > & \ 

£+l\ R ) 

c 

(54) 

B 0 s' - f K C / S ,r \ 

(k + nk 0 )R 0 A *t,s’j 0 d r r «<R 0 - r)x n>n J l+1 ( f- ) 


' (k + nk Q k 0 )^ - 


-- + A x 


n,n-l 


*Vl( P. 


E l+l,s r 


)• 


where usa has been made of w = (k + nk_)V, and (k + nk ) = k,./(l - V,/c). 

Ob o 0 b 

In the present analysis, it is assumed that the axial distribution 
function has the form 


G(C ) = * 


<C_ - y b "V b ) + 6 


2 ’ 


(55) 


where A is the characteristic spread in C about the mean value C « y, mV, . 

* z b b 
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We further assume that the characteristic spread A is small in comparison 
with Y b mV b . Substituting Eqs. (5) and (55) into Eqs. (28) and (41), 
we obtain the approximate expression 


“ 2 " + nk 0 Hk + n’k 0 )c 2 


n,n 


Y b ^o - U + nk Q )V b + i|k + nk 0 |i/Y b m] 2 
0 . Rq < r _< R c 


0 < r ' V 
(56) 


In obtaining Eq. (56), use has been made of Lc. (18). Making use of the 

definition of Budker's parameter in Eqs. (3) and (4), the term 4v/y,R 2 

b 0 

n 2 2 2 

in Eq. (56) can also be expressed as 4v/y, R“ « to ,/c , where w , « 

b U pb pb 

2 

4nnge /v^m is the plasma frequency-squared. 

Tlie condition for a nontrivial solution to Eqs. (50) and (51) 

is that the determinant of the coefficients c„ . and b be equal 

t,s i+l,s 

to zero. After some algebraic manipulation, we find that the TF. mode 
dispersion relation can be expressed as 


ui - (k + nkg)V b + i 


|k + nk 0 |a 


{ 4 - < k * "»o - V 2 - 


If |k + nk„ | A 

2 \ 

( “ - ( k * nk 0 )V b + 1 3 

- 4 VC ) 
3_2 / 

V Y b a 

Y b R 0 J 


4A 


q A 

2 vc 


YrR 2 
b c 


k Q (k + nk Q - kjj) 


VU 




P. 


-, E l a t+l,s R 0 
£ £ss'\ 


^tss-f^r 1 - 0 )!"- < k+nk 0 )V b +i 


|k + nk Q |A 


3 

Y b m 



-4^H '( 

3_2 Iss \ 

Y b R 0 


a £-tl ,S R 0 


where the coupling coefficient Q R , (a.., R n /K ) is defined by 

K.SS i|S v C 




*1+1.s 




'ft+l.s 


(£ + 1)2 j2 + lS + l,s> 


and the functions G is t(*) and H ffls ,(x) are defined by 


(57) 


(58) 
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2 J?(x) 

C £s' (x) ’ 7.J—3 


and 


(x" 

,2 


6 £.s‘> 


ll £ss' (x) 


W x) - J t (x)J t ^ 2 (x) 
G ts , ( x ) 


(59) 


(60) 


In Eq. (57), the subscript s and s' represent the radial mode numbers 
of the transverse and longitudinal perturbations, respectively. 

For small viggler amplitude (\ « 1), we investigate free electron 

laser stability properties for u and k + nkg near the simultaneous 

2 2 2 

zeros of the transverse dispersion relation, u. - (k + nkp - k^) c - 
2 2 2 

a i!.+l s C /R c = anc * t * le ^ on 8 iCu ^ ina 1 dispersion relation 


w - (k + nkp)V^ + i 


|k + nk Q |i 


y.m 


- 4 


v 3 R 2 

Y b R 0 


0 . 


(61) 


In this regard, making use of Eq. (61), the TE mode dispersion relation 
in Eq. (57) can be approximated by 


( H— - (k + nk - k ) 2 - ■■ ~ a 1 
|2 v 0 0 ; 2 I 

l c R J 

c 


k + nk 


it) — (k + nk^)V^ + i 




o^l ; 

) 


- 7. 


3 2 
Y b R 0 


44 


2 vc‘ 


Y.R 

b c 


k Q (k + nk 0 - k Q ) - 2 


0 E / ° H-1, s K Q 
^iss’V R 


In a similar manner, from Eqs. (53) and (54), we obtain the 


c 

(62) 


approximate TM mode dispersion relation, 

„2 1 


{: 


2 

2 ” 

- 4 


8 


(k + nkp - k Q ) 


£+l)S \ 


- (k + nk Q )V b + 1 


|k + nk Q |4 


12 


Y b m 


2 1 


3 2 

Vo 


4A 


2 vc 




k Q (k + nkjj - kp) 


"t+l.s 


M / 
Q £ss’ ( 


(63) 
B i+l,s K Q 


where the TK mode coupling coefficient Q, ,(B„R n /R ) I s defined by 
r 0 £ss £+l,s 0 c 

(64) 


<ss' (x) " G £s' (x)/J £ + 2 {6 £ + 1,s> ' 


and the function G fs ,(x) is defined in Eq. (59). 
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Figure 1 shows plots of C^ s ,(x) versus x obtained from Eq. (59) 

for (a) s ’ “ ®o i and (*>) gt “ Except in the case 

6. i = g n , , tlie plots of C , (x) for arbitrary 6 , are similar 

to those for g^ s' ” ^1 3 * n ^*6' ^(h). As shown in Fig. 1(b), the 

quantities cf . £ G. ,(x. ,) and x. , denote the maximum value of 
Is is is is 

G ts ,(*) and the corresponding value of x for a specified g^ ^,. For 
example, in Fig. 1, (x^,, G* gl ) = (0, 0.09) for B fL g , » g Q 1 and 
(x 0 , > G? i) “ (9.8, 0.064) for g , » g Shown in Fig. 2 are 

X- s £S £ , S 1 , J 

plots of (a) X £. 1 an< ^ 0>) tlie corresponding values of g| , = 

G , (x is , ) for several different values of the azimuthal and 

radial mode numbers i and s'. It is evident from Fig. 2(b) that 

C^ s , decreases rapidly with increasing values of the mode numbers i 

and s'. Moreover, we note from Fig. 2(a) that x {c _, can be approximated b 


X is' = 6 i.s- ' S ’ y 1 


(65) 


In this regard, for s' y 1, C* , can he approximated by 

K.S 

C L' - Ws-* - J i + i (fc i,s->- s ’ ' l - (66) 

r 

Shown in Fig. 3 are plots of (a) Q ,/G , for the TE mode 

£SS £K 

and (b) Q^ ss ,/G^ s , for the TM mode, obtained from Eqs. (58) and (64) 
respectively. Note that the curves in Fig. 3 are independent of tiie 
longitudinal radial mode number s'. Evidently, tlie values of Q ,/G , 

fSS £S 

and Q ,/G , increase with increasing values of azimuthal and 

£ss £s 

transverse radial mode numbers, l and s. After careful examination 
of Eqs. (58) and (64), we find that the maximum coupling between the 


transverse and longitudinal modes occurs for a value of R^/K^ given by 
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is' 

1 

t+l,s 


, /g 

is' 

1 p i+l,s 


R 0 /R c ' \ " *■**■ (6?) 

Equation (67) is valid only when x^, _< g for the TE mode, and 

x is- - <f W for the ™ raode - For x is' > Vi.s (TE >- or x t8 - > 6 i + l.B 
(TM), the maximum coupling occurs for R Q /R c “ 1. The maximum coupling 

coefficients corresponding to Eq. (67) can be determined from Figs. 2(b) 

and 3. For example, for (i,s,s') - (3,2,1), we determine that the 

maximum coupling coefficient and the corresponding ratio Rq/R • are 

given by (Q E gs , r q/ r c ) “ (1.607, 0.625) for the TE mode, and (Q*' ggl , Ro/R c ) 

(1.83, 0.52) for the TM mode. 

It is instructive to examine the present results for perturbations 

with the lowest mode numbers, i.e., (i,s,s') « (0,1,1), particularly 

for a beam-filled waveguide with Rq/Rc. ~ l. In this limit, from Fig. 1(a), 

we obtain G„ ,(a ) *■ 0.4 for the TE mode, and G„ ,(B, ,) “ 0.045 mode. 

Is 1,1 is 1,1 

We therefore conclude that the TE mode polarisation is the most unstable. 
Multiplying G ,(u. ,) - 0.4 by Q E ,/G, - 4.2 in Fig. 3(a), the 

£S 1 « X JLSS fcS 

E 

coupling coefficient is given by Qq^ “ 1.7. Assuming zero axial 

2 2 

momentum spread (a *= 0) and short axial wavelengths (k^R^ >:> 
the TE mode dispersion relation in Eq. (62) can be approximated by 


. 2 \ 

\ * \ 

I ^2 + nk 0 ’ ] 

<[w “ (k + nkp)V^] 2 - 3 

V c ' 

l *b 1 


, 2 2 .2 
ul pb k 0 ’ 


( 68 ) 


for the U.s.s’) = (0,1,1) perturbation and R./R 

0 c 


1. Equation (68) 


is similar in form to the result obtained by Davidson and Uhm for 


a uniform density beam with infinite cross section. In particular, 
the constant numerical factor on the right-hand side of Eq. (68) is 
equal to 3.4, whereas in Ref. 3 the constant numerical factor is equal to 8. 
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Finally, we have investigated detailed stability properties 
by solving the dispersion relations in Eqs. (62) and (63) numerically 
for a broad range of system parameters. Defining the normalized 
Doppler-shifted eigenfrequency by 

ii - [w - (k + nk 0 )V fc )/k 0 c , (69) 

we calculate the normalized growth rate 0^ = lmO from Eqs. (62) and 

(63). Shown in Fig. 4 are plots of the normalized growth rate 

versus (k + nk Q )/k 0 for U,s,s') = (3,2,1), k^ = 10, y b * 10, v/y b “ 

2 

0.02, and A • 0.01, with (a) R../R = x,,/u. . for the TE mode, and 

U c jl 4,2 

(b) R n /K « x.,/6, - for the TM mode. For these optimized choices of 
U c 31 4,2 

R 0^ R c’ t * 1e i nsta bili c y growth rate for the TM mode is comparable to 

that for the TF. mode. Moreover, the growth rate is reduced substantially 

by introducing a small amount of axial momentum spread (A/y b mc £ 0.01). 

We conclude this section by pointing out two areas in which the 
analysis can be extended. First, the restriction to very short ;wave- 
length perturbations [Eq. (46)] can be removed in a relatively straight¬ 
forward manner. Second, paralleling the self-consistent theoretical 

g 

formalism developed in previous studies, the stability analysis can 
also be carried out without making the approximation that the transverse 
perturbations are represented by the vacuum waveguide fields. 
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V. CONCLESIONS 

In this paper, we have examined the free electron laser instability 

for a solid relativistic electron beam propagating in the combined 

transverse wiggler and uniform axial guide field given in Eo. (1). 

The analysis was carried out within the framework of the linearized 

Vlasov-Maxwell equations. The equilibrium (Sec. II) and stability (Secs. 

Ill and IV) properties were investigated in detail for the choice of 

distribution function in which all electrons have the same value of 

the linear combination of transverse and helical invariants, C, - 2v,nu), C, , 

b b h 

and a Lorentzian distribution in the axial invariant [Eqs. (5) and 

(18)]. One of the most important conclusions of this analysis is that 

the maximum instability growth rate for a solid electron beam is comparable 

to that of i hollow beam win sin-ilai parameters. horeover, it is also 

found that the maximum growth rate occurs at a value of Rq/R^ corresponding 

to Rn/R_ * x , /a.,, for TE mode perturbations, and F.,,/R « x„ 

^ c £s t+l,s 0 c Js J+l,s 

for Tfl mode perturbations. For these optimized values of R./R , the 

U c 

instability growth rate for the TK inode is comparable to that for the TL 
mode. Moreover, the growth rate is substantially reduced by introducing 
a small amount of axial momentum spread (i/y^mc ^ 0.01). 
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FIGURE 2 PLOTS OF (a) x gs - AND <b) THE CORRESPONDING G$ s > = G gs - (x gs <) FOR SEVERAL VALUES 
OF AZIMUTHAL AND RAIDAL MODE NUMBERS, i. AND S’. 
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(a) TE mode 
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(b) TM mode 
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A 


FIGURE 3 PLOTS OF (a) Q| s */G gs - AND (b) Q^ s '/G es ' 

EQS. (58) AND (64) FOR SEVERAL VALUES OF 8 AND S. 
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FIGURE 4(a) PLOTS OF NORMALIZED GRWTH RATE VERSUS (k + nk 0 )/ko FOR (e,8.s’> - (3,2,1), 7 b “ 10, 
vly b - 0.02, AND A 2 = 0.01, WITH (a) Rq/R,. - x 31 /c* 4 2 F0R THE TE M0DE - 






FIGURE 4(b) R 0 /R c = x 31 /0 4 2 FOR THE TM MODE 
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APPENDIX A 

LONGITUDINAL PERTURBATIONS FOR THE FREE ELECTRON LASER INSTABILITY 


In this Appendix, we investigate properties of the longitudinal 
perturbations about an electron beam propagating through a cylindrical 
waveguide with radius R^. In the present analysis, it is assumed that 
the perturbations have short wavelength with 

U*1 - |<k + nk Q ) 2 - w 2 /c 2 | » 1 /Rq , (A.l) 

which can also be expressed as 

q n R 0 " (1 + V'^b^O ~ 1 * (A ‘ 2) 


for the frequencies of interest for free electron laser applications. 
Equation (A.2) is easily satisfied in parameter regimes of present 
experimental interest. In the limit of a small uiggler amplitude (A -* 0), 
we obtain the longitudinal eigenvalue equation. 


/ 1 3 3_ l 2 \ ‘ (n) . , 

( r 3r r 3r “ 2 ~ q n ) E z,i (r) 




3_,2 ’ 


[w - (k + nk 0 )V fa + ijk + nk 0 |d/i b m] 
from Eqs. (33), (45), (46), and (56). In Eq. (A.3), C(x) 


step function defined in Eq. (52), and w 


pb 


4vc 2 /Y b Ko is 


frequency-squared. 

For notational simplicity, we define 


(A.3) 


is the Heaviside 
the plasma- 


i/(r) = E<"j(r) . 

Z | X 


(A. 4) 


Inside the electron beam (0 < r < Rq) > Eq- (A. 3) can be expressed as' 


12 
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{$h r h --2 + r2 ) 6 * t(l) ■ 0 • °-‘ r <R o • 


(A.5) 


where 


2 . 2 
T = q 


2 / 2 
M pb /Y b 


3_,2 


" [[“ - (k + nk 0 )V b + i|k + nk 0 |A/v b m] 


1>. (A.6) 


Outside the electron beam (R. < r < R ), Eq- (A.3) reduces to the free- 

u c 

space eigenvalue equation 

( rS7 r IT-4 - *'n )4 * t(r) * ° ■ R 0 ‘ r <R c ; (A ‘ 7 > 


The solution to Eq. (A.5) that remains finite at r ” 0 is 


6^ n (r) - i t J t (TR) , 0 r < R Q , (A. 8) 

where J^(x) is the Bessel function of the first kind of order £, and 

* 2 2 

$ is a constant. Noting a R' >> 1 in Eq. (A. 2), »’e can express the 
£ 'n U 

solution to Eq. (A.7) as 


6^ ut (r) - C[l £ (q n r)K £ (q n R c ) - K (q^) I £ (q^) ], ^ < r < » C ,(A.9) 

where 1^ and are modified Bessel functions of order i , and C is a 
constant. 

The boundary conditions on iif^r) at the surface of the electron 
beam are given by 

■ K«w 0 • <*•“> 

and 

l £) ' Sr)4 »‘in 3 r.R 0 ■ * (A ’ U) 

from Eq. (A.3). Substituting Eqs. (A.8) and (A.9) into Eqs. (A.10) and 
(A.11) gives 


A-2 
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TR 


j ;< tr 0 > 

0 J £ (tr 0 ) 


h( % ) 


(A.12) 

- p 1 i (, .W , » R c ) ' Wc>Wo> 

' V ° V'nW^V ’ V q n P c )K £ (q n R 0> ' 
where the "prime" notation denotes derivative with respect to the complete 
argument of the Bessel function, e.g., J'(TR 0 ) - [dj (x)/dx) . The 

£ U £ X~TRq 

expression for the longitudinal wave admittance h(q n ) in Eq. (A.12) 

can be simplified in several limiting cases, including short wavelength 
2 2 

perturbations with |<J n RQ | >> 1. In this case, h(q^) can be approximated by 


h(q n ) = " q n R 0 COthq n (R c " V ’ (A ' 13) 

and Eq. (A.12) reduces to 
J 1 (TR ) 

- TR 0 TTfR^j" “ q n R 0 COthq n (R c " R 0 ) ' (A.1A) 

Evidently, the right-hand side of Eq.(A.14)is a very large positive 
number, and the lovest-order longitudinal dispersion relation (for t -> 0) 
can be approximated by 


V TR 0 ) -■ 0 ’ 


(A.15) 


where T is defined in Eq. (A.6). It follows from Eq. (A.15) that 


2 2 2 

T R 0 “ 6 t, s ’ * s ’ " l ’ 2 . 


(A. 16) 


where fi , is the s’th aero of J (x) «= 0. In this regard, Eqs. (A. 8) 
£ * S £ 

and (A.9) can be approximated by 




> £,s ,J £ (B £,s ,r/ V ’ ° - r * R 0 ’ 


otherwise , 


(A.17) 


where $ , is a constant. Substituting Eq. (A.6) into Eq. (A.16) 

£* 8 


A-3 
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and making use of Eq. (A.2), we obtain the longitudinal dispersion relation, 


2 2 


u - (k + nk„)V, + i 
U b 


|lc + nkfili'l w , 

-_ _E^ . o 

3 2 

*h n 


(A.18) 


where the term proportional to 6* g , has been neglected. 


A-4 
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